Gravity-induced resonances in a rotating trap 
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It is shown that in an anisotropic harmonic trap that rotates with the properly chosen rotation 
rate, the force of gravity leads to a resonant behavior. Full analysis of the dynamics in an anisotropic, 
rotating trap in 3D is presented and several regions of stability are identified. On resonance, the 
oscillation amplitude of a single particle, or of the center of mass of a many-particle system (for 
example, BEG), grows linearly with time and all particles are expelled from the trap. The resonances 
can only occur when the rotation axis is tilted away from the vertical position. The positions of the 
resonances (there are always two of them) do not depend on the mass but only on the characteristic 
frequencies of the trap and on the direction of the angular velocity of rotation. 

PACS numbers: 03.75.Kk, 45.30.+S, 45.50.-j 



I. INTRODUCTION 

The main purpose of this work is to expose the role of 
gravity in the dynamics of particles in a rotating trap. 
In addition, we present a complete analysis of the sta- 
bility regions for a rotating trap in 3D. We prove that 
in the generic case there are three separate regions of 
stability with different characteristics. Gravity induced 
resonances are relevant only if they occur in the regions 
of stability, otherwise, they are swamped by the expo- 
nential behavior of trajectories. 

Harmonic traps are often used in optics and atomic 
physics (especially in the form of TOP traps P| in the 
study of Bose-Einstein condensates) and yet a complete 
theory of these devices has not been developed. The so- 
lution to the problem of a trap rotating around one of the 
trap axes is effectively two-dimensional and its solution 
has been known for at least hundred years. In the clas- 
sic textbook on analytical dynamics by Whittaker Q we 
find a solution of a mathematically equivalent problem of 
small oscillations of "a heavy particle about its position of 
equilibrium at the lowest point of a surface which is rotat- 
ing with constant angular velocity about a vertical axis 
through the point" . The quantum-mechanical counter- 
part of the Whittaker problem has also been completely 
solved 0, and the statistical mechanics of a classical 
gas was studied in |^. 

In this work we present a complete solution to the 
problem of the motion of a particle moving in a most gen- 
eral anisotropic rotating harmonic trap in 3D and in the 
presence of gravity. This is an exactly soluble problem 
but technical difhculties apparently served so far as a de- 
terrent in developing a full description. A full description 
of the particle dynamics in a rotating anisotropic trap in 
three dimensions so far has not been given despite a new 
significance of this problem brought about by experimen- 
tal and theoretical studies of Bose-Einstein condensates 
and the accompanying thermal clouds in rotating traps 
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lid llU Il2j . Explicit formulas describing the 
complete mode structure in the three-dimensional case 
are indeed quite cumbersome because we deal here 
with third-order polynomials and on top of that they 
have rather complicated coefficients. However, many im- 
portant features may be exhibited without straining the 
reader's patience. In particular, we can identify various 
stability regions for an arbitrary orientation of the angu- 
lar velocity and we can give conditions for a resonance. 

The standard arrangement |lllllllli|l3,|ll|llis 

to choose a vertical axis of rotation of the trap. Slight 
tiltings of this axis were introduced to excite the scissors 
modes 0,0,0]. However, for such very small tilting 
angles the effects described in the present paper would 
not be noticeable. In the case of a vertical axis of rota- 
tion, there are no resonances. The only effect of gravity 
is a displacement of the equilibrium position. The situa- 
tion completely changes and new phenomena will occur 
when the axis of rotation is tilted away from the verti- 
cal position. In this generic case, for every anisotropic 
three-dimensional trap there exist two (not three as one 
might expect) characteristic frequencies at which reso- 
nances occur. The motion in a trap that is rotating at a 
resonant frequency will become unbounded and all par- 
ticles will be expelled from the trap. The position of the 
resonance does not depend on the mass but only on the 
characteristic frequencies of the trap and on the direction 
of the angular velocity. 

All our results are valid not only for a single par- 
ticle but also for the center of mass motion in many- 
body (classical or quantum) theory since for all quadratic 
Hamiltonians the center of mass motion completely sep- 
arates from the internal motion [T^ IT^ [l8l|. Therefore, 
a trap rotating at the resonant frequency will not hold 
the Bose-Einstein condensate. Owing to the linearity of 
the equations of motion for a harmonic trap, all conclu- 
sions hold both in classical and in quantum theory. A 
resonant behavior caused by an application of a static 
force may seem counterintuitive, but it is explained by 
the fact that in a rotating frame the force of gravity acts 
as a periodically changing external force. 
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II. EQUATIONS OF MOTION 

The best way to analyze the behavior of particles in a 
uniformly rotating trap is to first perform the transfor- 
mation to the rotating frame. In this frame the harmonic 
trap potential is frozen but the force of gravity is rotat- 
ing with the angular velocity of the trap rotation. In the 
rotating frame the Hamiltonian has the form 



zm 



rflp + —rV-r 



mr-g{t). 



(1) 



The potential matrix V is symmetric and positive def- 
inite. The eigenvalues of this matrix are the squared 
frequencies of the oscillations in the non-rotating trap. 
The angular velocity matrix f2 is related to the compo- 
nents of the angular velocity vector through the formula 
^ik = £ijk^j- The vector of the gravitational accelera- 
tion g{t), as seen in the rotating frame, can be expressed 
in the form 



g{t) = g\\ + g± cos{nt) - (n x g^) s'm{nt), 



(2) 



where n denotes the direction and J7 denotes the length 
of the angular velocity vector ft. The parallel and the 
transverse components of the gravitational acceleration 
vector g = g{0) are defined as g|| = n{n g) and g±^ = 
g — n{n-g), respectively. Note that the time-dependent 
part vanishes when the rotation axis is vertical. 

The equations of motion determined by the Hamilto- 
nian have the following form 



drjt) ^ p{t)_ 

dt m 
dp{t) 



n-r{t), 



(3a) 



dt 



^mV-r{t) -h-p{t) +mg{t). (3b) 



These equations describe an oscillator in a rotating frame 
displaced by a constant force (the longitudinal part of g) 
and driven by a periodic force (the transverse part of g). 
It is convenient to rewrite the expression (|2Jl as a real 
part of a complex function 



g(<)-5R(sf||-f (Sfi+z(nxg^))e'"*) 
In compact notation Eqs. (|3a|l have the form 
dn{t) 



dt 

where 

n{t) 

G\\ = m 



Min)-nit) + ^{gn+g^e 



frit) 



9\\ 



, g± = m 



-mV -Cl 


g± +i{n X gj_) 



(4) 
(5) 

(6) 
•(7) 



We shall now replace the equations of motion by their 
complex counterpart 



dW{t) 
dt 



Min)-w{t) + g»+g_Le 



int 



(8) 



The physical trajectory in phase space is described by the 
real part of the complex vector yV(i). Let us introduce a 
basis of six eigenvectors of M{fl) 

Min)Xk^iuJkmXk, fc = l,...,6 (9) 
and expand W(t) and g{t) in this base as follows 
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Wit) = ^a'^(i)e*'^'=(")* A-fe, (10) 

6 6 

^11 = T.^\\^k, g^ = Y.^iXk. (11) 



fe=l 



fc=l 



Owing to a simple block structure of A4{fl), the basis 
vectors Xk can be determined by reducing effectively the 
problem to three dimensions. We use this method in the 
Appendix B to determine the resonant solution. 

The equation of motion (|SJ) can be rewritten now as a 
set of equations for the coefficient functions a'' (t) 



da^{t) 
dt 



= 1,...,6. 

(12) 

It is clear now that the mode amplitude ak{i) will grow 
linearly in time — the signature of a resonance — when- 
ever either one of the two terms on the right hand side 
becomes time independent. This happens to the first 
term if one of the frequencies u;fe(f2) vanishes but the 
corresponding coefficient does not vanish. This case 
is not interesting, since it means that we are just at the 
border of the lower instability region and the trap is not 
holding particles, as discussed in the next section. The 
second term becomes time independent when the angular 
velocity of trap rotation satisfies the resonance condi- 
tion 17 = ijJk{^) and, of course, 7^ 7^ 0. This resonance 
is different from a resonance in a standard periodically 
driven oscillator. In the present case the characteristic 
frequencies of the trap depend on the frequency O of the 
driving force. Therefore, the position of the resonance 
has to be determined selfconsistently. A full description 
of these gravity induced resonances requires the knowl- 
edge of the behavior of a;fe(Sl)'s as functions of Vl. In 
particular, it is important to know whether a resonance 
occurs in a region where the system undergoes stable os- 
cillations. This will be discussed in the next Section. 



III. REGIONS OF STABILITY 

The stability of motion for a harmonic oscillator is 
determined by the values of its characteristic frequen- 
cies w — the roots of the characteristic polynomial. In 
the present case, these frequencies are determined by the 
characteristic equation for the matrix M.{^) 



Det 



^M{n)~iuj^ 



= 0. 



(13) 
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The characteristic polynomial is tri-quadratic 



(14) 



where the coefhcients A, B i C can be expressed in a 
rotationally invariant form jT3| 



A = -2n^ - Ti-{V}, 

B = Q.^ + {f{in-V-n~TT:{V})- 



Tr{V^}2-Tr{y2} 



C = Q?-{i:Y{V}-Q?)n-V-n~^fn-V^-n~V>ei{V}. (15) 

Stable oscillations take place when all characteristic tj's 
are real. This means that all three roots of the polyno- 
mial Qix) must be real and positive. Without rotation, 
when 17 = 0, the three roots of Q{x) are equal to the 
eigenvalues of the potential matrix V . We have then a 
simple system of three harmonic oscillators vibrating in- 
dependently along the principal directions of the trap. 
As n increases, our system will, in general, go through 
two regions of instability: the lower region when one of 
the roots of Q{x) is negative and the upper region when 
two roots are complex. We shall exhibit this behavior by 
plotting the zero contour lines of Q{x) in the riX'Pian^- 
We assume that the trap potential and the direction of 
rotation are fixed and we treat the characteristic polyno- 
mial Q{x) as a function of and x only- Contour lines 
representing the zeroes of Q{x) i^i the generic case are 
shown in Fig.^ There is a region of f2, where only one 
real root of Q{x) exists. However, this region is bounded, 
so for sufficiently large the system is always stable. 

It has been argued in Ref. that there is always 
a region of instability when one of the roots of Q{x) is 
negative. The corresponding modes grow exponentially 
with time. As seen in Fig. this region of instability 
is bounded by the two values ^1,2 at which the curve 
crosses the vertical axis. These values are given by the 
zeroes of C, treated as a biquadratic expressions in f2 



2^ ' 



(16) 



where a = n V n, b — Ti{V}n-V -n — n-V^-n, and c = 
Det{y}. Since a, 6, and c are positive and iP' > Aac, both 
values r^i and O2 are real. A degenerate case is possible, 
when Jli = then the region of instability shrinks to 
zero. In order to determine, when this can happen, we 
may use the (explicitly non-negative) representation of 
the discriminant 6^ — 4ac given in Ref. . Assuming for 
definitness that 14 < < T4, we find that this happens 
in two cases: when the trap is not fully anisotropic {Vx — 
Vy or Vy = Vz) ot: the axis of rotation lies in the a;z-plane 
and its azimuthal angle satisfies the condition sin^ 9 = 
(1 — Vx/Vy)/{1~ Vx iVr. ) ■ The second possibility has not 
been noticed in Ref. [iq. 

Graphical representation of the solutions for Ux = 1, 
Uy = Uz = is shown in Fig. |21 In this plot only one in- 
stability region is present, where one of the roots of 114|) 



is negative. Owing to the stabilizing effect of the Coriolis 
force, for fast rotations the system becomes again stable. 
We can also see that one of the characteristic frequencies 
remains constant, it does not vary with Q. It is so, be- 
cause the rotation does not influence the motion in the 
direction of the rotation axis. This is the degenerate case 
described by Whittaker and thoroughly studied in con- 
nection with the BEC traps 3, 4, 5, lOj. In this case 
there is only one instability region, where the square of 
the frequency is negative. When the direction of angu- 
lar velocity is not parallel to one of the axes of the trap, 
a second kind of instability appears. In addition to the 
region of (purely exponential) instability, described be- 
fore, when one root of Q{x) was negative there is also, in 
general, an additional region of (oscillatory) instability 
where two roots are complex. Since the coefficients of 
the characteristic polynomial are real, the square of the 
second frequency is complex conjugate to the first one. 
In this case the instability has the form of expanding 
oscillations. In Fig. |3| we show how the second kind of 
instability develops when the rotation axis is tilted away 
a little bit from the z direction. 

We will show now, that the second kind of instabil- 
ity does not occur in the most common case, when the 
rotation vector is along one of the axes of the trap. In 
particular, it can never happen in a two dimensional trap. 
Without any loss of generality, we may assume that the 
direction of rotation is along the z axis. In this degen- 
erate case the characteristic polynomial (|14|) factorizes 
because the motion in the z direction is not influenced 
by rotation, 

Q{x) = {x~Vz) 

x{x^ - x{2n^ + Vx + Vy) + n^(Vx + Vy) + VxVy). 

(17) 

The factor quadratic in x has real zeros if its discriminant 
A is positive and this is, indeed, the case since 

^ = %^l'{Vx+Vy) + {Vx-Vyf >Q. (18) 

It follows, however, from the topology of the curves rep- 
resenting the zeros that the second kind of instability al- 
ways exists if the direction of rotation does not coincide 
with one of the axes of the harmonic trap. 



IV. GRAVITY INDUCED RESONANCES 

There are, in general, two resonant values of 17. They 
occur when x = ^^ ■ In this case, from Eq. H14|) we get a 
biquadratic equation for the resonant values of f2: 

DVf^ + £^1^ +F ^Q, (19) 



where 

D = -2{Tt:{V} - n-V -n) 
Tr{y}2-Tr{l/2} 



E = 



F = -Det{y}. 



TY{V}n-V -n- n V'^ n, 



(20) 
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FIG. 1: The contour hues in this figure represent the zeroes 
of the characteristic polynomial 1141 for Vx — 1, Vy — 2, 
Vz — 3 and Ux = l/\/3, riy — Uz = l/y/S plotted 

as functions of the magnitude of angular velocity Q and the 
square of the characteristic frequency x = t^^- In addition to 
the lower region of instability (between dashed lines) where 
one of the roots of 1141 1 is negative, there is also an upper 
region (between dotted lines) where there exists only one real 
root (determined by the continuation of the line that begins 
at X = 3) not seen in this figure because it corresponds to 
a very large value of x)- Since there is no absolute scale of 
frequencies involved in the analysis of stability regions, we 
have chosen the lowest trap frequency as a unit. 
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FIG. 2: This plot is for the same trap as in Fig. Q but for 
the rotation around the trap axis Ux = 0, Uy = 0, Uz — 1. In 
this degenerate case there exists only one region of instability, 
when one of the roots of 11411 is negative. 

The properties of the two solutions j_ of Eq. H19|l are de- 
scribed in the Appendix A. We may also analyze the reso- 
nance condition graphically by superposing the parabola 
X = (shown as a thick line in Figs. 2K7|) on the plots 
of characteristic frequencies. The lower value is more in- 
teresting because, as shown in the Appendix A, it always 
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FIG. 3; This plot shows how the degenerate case (Fig. I^J 
merges with the general case (Fig. 0. Here we have chosen 
the following values of the parameters: Vx = 1, Vy = 2, Vz ~ S 
and Hx = sin(l/10), Uy = 0, riz = cos(l/10). 




FIG. 4: Vx = I, Vy = 2, Vz = i Ux = sin(^), Uy = 0, 
Uz — cos(^). Two horizontal lines enclose the lower region 
of instability. Both resonant frequencies lie in the lower region 
of stability. 

falls into the range of the first region of stable oscillations 
where we would expect the confinement of particles in the 
trap. The upper value may fall in the lower region of sta- 
bility (Fig.0J), in the lower region of instability (Fig.EJ, 
or in the higher region of stability (Fig.EJ. A degenerate 
case is also possible (cf. Fig. [Tl when the two resonance 
frequencies coincide but this happens only under very 
special circumstances (see Appendix A). The existence 
of only two resonant frequencies (we would expect three 
resonant frequencies for a three dimensional oscillator) 
clearly shows that we are dealing here with a more com- 
plex dynamical system than a simple driven harmonic 
oscillator. 

The essential difference between the resonant and the 
nonresonant behavior is illustrated in Figs. 181101 The role 
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FIG. 5: K = 1, Vy = 2, V, = Z = sin(f), Uy = 0, 
= cos(-|). Higher resonant frequency lies in the lower 
region of instability. 



FIG. 7: In the degenerate case, when two resonant frequen- 
cies coincide, their common value is also equal to the lowest 
frequency of the trap. In this case the trap parameters are: 
14 ^l,Vy = 10/3, 14 = 5, and = 1. 




1 

FIG. 6: 14 = 1, 14 = 2, 14 = 3 = sin(^), Uy = 0, 
Tiz = cos(^). Resonant frequencies lie in two different regions 
of stability. 



of gravity is best seen by comparing Fig. |H1 (with grav- 
ity) and Fig. ^2 (gravity switched off) . The calculations 
were performed in the coordinate system in which the 
potential matrix V is diagonal and it was assumed that 
the force of gravity at t = is directed along the z axis 
of the trap. All four figures were generated for the same 
trap and under the same initial conditions: the particle is 
initially placed in the center of the trap and it is given the 
initial velocity of 1 cm/s in the x direction. These simple 
initial conditions result in the excitation of all the modes 
of the oscillator. However, at resonance the mode grow- 
ing linearly with t dominates the time evolution. This 
can be seen by comparing, for the same parameters of 
the trap and at the resonant frequency ri/27r = 6.49421 
Hz, the motion depicted in Fig. |S1 that was calculated 
numerically, with the motion depicted in Fig. 1121 The 
second plot was obtained from an analytic solution of the 



equations of motion JHJ), that contains only the resonant 
part. This analytic solution is derived in the Appendix 
B. Both plots are essentially the same except for very 
small differences that are due to the fact that the ini- 
tial conditions assumed in Fig. [S] require some admixture 
of the nonresonant solutions while the analytic solution 
does not contain any nonresonant pieces. 

Gravity induced resonances are significant only if they 
occur in the region of stability. In all regions of insta- 
bility, where the trajectories exhibit exponential growth, 
the resonances cannot be detected. 



V. CONCLUSIONS 

We have analyzed the stability of motion in an 
anisotropic, rotating harmonic trap in 3D. We have found 
that, in general, there are three regions of stability. The 
second region and the third region merge only when the 
rotation axis coincides with the one of the trap axes. We 
have demonstrated the presence of resonances in a rotat- 
ing harmonic trap subjected to the force of gravity. The 
nature of these resonances is different from the standard 
forced harmonic oscillator since the resonant frequencies 
must be determined selfconsistently. The resonances oc- 
cur at two rotation rates but they exist only when the 
rotation axis is tilted away from the vertical direction. 
However, the rotation axis can still be directed along one 
of the trap axes. The lower resonance always falls in 
the region of stable oscillations where as the higher reso- 
nance may fall in the lower region of stability, in the first 
region of instability, or in the higher region of stability. 
Resonant rotation rates depend solely on the properties 
of the trap and not on the masses of particles. The res- 
onances cause the escape of a particle, or the center of 
mass for a collection of interacting particles, from the 



FIG. 8: The trajectory of a particle in the coordinate frame 
rotating with the trap. The characteristic frequencies of the 
trap in the x, y, and z directions are 10 Hz, 15 Hz, and 20 
Hz, respectively. The angular velocity vector has the direc- 
tion (1, 1, 1) and its length VL/2-n = 6.49421 Hz satisfies the 
resonance condition, as calculated in the Appendix A. The 
distances in this plot are measured in centimeters. The direc- 
tion of the gravitational force {g — 9.81ms~^) is assumed to 
coincide at t = with the z-axis of the trap. We can see that 
the amplitude of oscillations increases linearly with time. 




FIG. 9: The trajectory of a particle under the same conditions 
as in Fig. |S1 except that the value of the angular velocity 
Q./2n = 4 Hz is below the resonance. The scale in this plot is 
reduced by a factor of 10 as compared to Fig. |S1 



FIG. 10: The trajectory of a particle under the same condi- 
tions as in Fig.|5] except that the value of the angular velocity 
r2/27r = 9 Hz is above the resonance. The scale in this plot is 
reduced by a factor of 10 as compared to Fig. |S1 




FIG. 11: The trajectory of a particle under the same condi- 
tions as in Fig. |S1 except that the gravitational field has been 
turned off. The angular velocity Q,/2n — 6.49421 Hz has the 
resonant value. The scale in this plot is reduced by a factor 
of 100 as compared to Fig. |H] Clearly, there is no sign of any 
resonant behavior. 



trap. In view of recent experiments in which the axis 
of rotation was tilted away from the trap axis [T3. ITU , 
it should be possible to confirm experimentally the exis- 
tence of the gravity-induced resonances for Bose-Einstein 
condensates in rotating traps. In particular, the presence 
of resonances makes it possible to quickly expel all par- 
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tides from the trap. 
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Ministry of Scientific Research and Information Technol- 
ogy under the Grant No. PBZ-MIN-008/P03/2003. 



APPENDIX A: CALCULATION OF RESONANT 
ANGULAR VELOCITIES 



FIG. 12: The trajectory of a particle described by the analytic 
formulas given in the Appendix B. The trap parameters and 
the angular velocity are the same as in Fig. |H| 



J 



The resonant values of the angular velocity are the 
ones that coincide with a characteristic frequency of the 
trap. Thus, the resonant values of O are obtained by 
solving the biquadratic equation (|19|l 



Q{n^) = -2nHTY{V} - n-V-n) + f]M Tr{T^ } , ^^^^y^ ^ ^ ^ ^y^ 



TriVjn-V-n-n-V'-n - Bet{V} = 0. (Al) 



In the coordinate system aligned with the principal axes 
of the trap, the solutions of this equation have the form 



J 



ni 



(1 - nl/2)VyV, + • • • ± - nl/2)VyV, + ...)■'- 2V,VyV4{l - nl)V, + ...) 



2(1 



)X4 



r 



(A2) 



where Vx,Vy,Vz denote the diagonal elements of V and 
the dots stand everywhere for two additional terms ob- 
tained by the cyclic substitutions x y z x. These 
solutions are always real since the discriminant A — the 



J 



expression under the square root — is never negative. 
This is seen from the following representation of A as a 
sum of two nonnegative terms (for definitness, we have 
assumed here that Vx < Vy < Vz) 



A = {il-nl/2)VyV,-{l + nl/2)VxVz-{l+nl/2)VxVyy 



Wx {nliVz - VxKVyVz + V^/2) + nliVy 

I 



Vx){VyVz + V,y2)) 



(A3) 



This representation of A can also be used to find the into one. The second term vanishes (for a nondegener- 
necessary condition for the two resonant values to merge 
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ate trap) only if = ±1 and the first term then van- The lower value ri_ of the resonant angular velocity never 
ishes when the trap frequencies satisfy the condition (cf. exceeds the lowest critical angular velocity at which the 
Fig. 13 system becomes unstable. To prove this, we may com- 

pare the expression lAip with the following expression 
^ — J_ _l_ J_ determines the critical value | 16j 



Q(0) = -n'^ n-V-n + n^{TT{V} n-V-n - n-V^ -n) - Bct{V} = 0. (A5) 

I 



Both these expressions may be represented by inverted 
parabolas that cross the y-axis at the same negative value 
— Det{y}. Since the derivative of (|A1|) with respect to 
fi^ at the crossing point is larger than the derivative of 
(|A5I) . the parabola (|Aip is steeper at the crossing point. 
Therefore, it must cross the cc-axis at a lower value of fl^ 
and the lower resonant value lies below the boundary of 
the stability region. 

APPENDIX B: ANALYTIC RESONANT 
SOLUTION 

A resonant solution is the one that has in the mode 
expansion H1U|I only the term oscillating with the reso- 
nant frequency Q. The amplitude of these oscillations 
is growing linearly in time. In order to find an explicit 
form of this solution, we substitute into the equations of 
motion JSJl the following Ansatz 

Wr = iAt + B)e'^* +C, (Bl) 

where B, and C are six-dimensional vectors. Upon 
substituting this form into the equations of motion and 
comparing the terms that have the same time depen- 
dence, we obtain the following equations for the three 
unknown vectors 

^M{n)-A = iflA, (B2) 

'-M{n)-B = inB + A-gA_, (bs) 
^M{n)-c = -g\\. (B4) 



The first equation says that the vector A is directed along 
the resonant eigenmode. From the second equation we 
can determine the length of this vector. The vector B 
can be determined from the second equation only up to 
a component along A. This is so, because by changing 
the origin of the time scale, we may always add such 
a component to B. Finally, the vector C is determined 
from the third equation. It is now a matter of pure al- 
gebra to find these three vectors. Six-dimensional prob- 
lems are quite cumbersome but we may replace them 
here by their three-dimensional counterparts, owing to a 
simple block structure of the matrix 'M(fJ). In addition, 
let us note that we need only the upper three compo- 
nents of the solution Wr(i) to determine the trajectory. 
We shall call them a, 6, and c, respectively. Thus the 
physical trajectory for the resonant solution will have the 
form r{t) — di{{at + b)e*^^* + c). Eliminating the lower 
three components, we obtain the following set of three- 
dimensional equations 



N{n)-a^O, (B5a) 
N{fl)b = iflb + a-h, (B5b) 
7V(0)-c = -g||, (B5c) 



where iV(ti>) = uj^ —2iu;d—d^ — V and h = g±+i{nxg^). 



(up' + — n^) — Vx —Vi^nxTiy — 2iijjVtnz) —Vl^nxTiz + 2iijjVlny) \ 
-n'^n^ny + 2iuj^nz) lu^ + fl'^ {1 - nl) ~ Vy -Q'^nyTiz ~ 2iujQnx) \ . (B6) 
—VPnxnz — 2iuj^lny) —fl'^nyUz + 2iujnnx) lo'^ + fl"^ {1 — nl) — Vz / 



The best way to find solutions of Eqs. (|B5(I is to expand matrix N{ft) 
the vectors a and b in the basis of the eigenvectors of the 

a = aQBQ + a+e^ + a^e^, (B7a) 

b = /3oeo+/3+e++/3_e_. (B7b) 
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Upon substituting these formulas into Eq. ljB5a|l one finds 
that a_|_ and a_ vanish, Pq is arbitrary, and the remain- 
ing three coefficients can be calculated from Eq. l|B5b(l . 
Everywhere in this Appendix f2 stands for one of the two 



resonant frequencies given by the expression IIA2|l . The 
characteristic polynomial of the matrix N{Q) is 



+ A2(K + Vy + V,- 5r!2) + 4A(r!4 _ sn^^V, + Vy + K) + VyV, + + V^Vy). (B8) 

I 

One of the eigenvalues vanishes and the remaining two 
eigenvalues A± have the form 



A± = n =•= 



V.' - 2VyV, 



(B9) 



where the dots have the same meaning as in the Ap- 
pendix A. Knowing the eigenvalues, we may write down 
explicit expressions for the projectors Po,P+, and P+ 
corresponding to the eigenvectors of N{n), 



Po 

P+ 



{N{n) - x+)iN{n) - A_) 

A+A_ ' 

N{n){N{n) - A_) 

A+(A+-A_) ' 

N{n)iN{n) - A+) 

A-(A--A+) ■ 



(BlOa) 
(BlOb) 
(BlOc) 



The eigenvectors bq and e± may be obtained by acting 
with the projectors on any generic vector (not an eigen- 
vector of M(f2)). We shall choose h = g± + i{n x g±) 
as this vector because h appears already in Eq. IjBScp . 
Thus, the (unnormalized) eigenvectors, needed in the ex- 
pansions (|B7a|) . can be written in the form 



eo = Pq-H, 
e+ = P+h, 
e = P^h. 



(Blla) 
(Bllb) 
(Bllc) 



Taking into account the orthogonality of the eigenvectors, 
we may find the coefficients and j3± that lead to the 
following explicit expression for the vectors a, b 



a = 



eo, 



b = 



»(eo'go) (e+-^-eo) 
{el-{in + Cl)-eo) [e^-e 

i(ej- eo) (e'L-Cl-eo) 



{el-{ifl + n)-eo) 



(B12a) 



-1 



-.(B12b) 



Finally, the vector c is found by solving the equations 
M{Q)-c — —n{n g) and it has the following components 



V^VyV, 



^^{nlV,{Vy + V.\ 



nlVyiV, + V,)+nlV,{V, 
I 



Vy)) + n^{nlV,+nlVy + nlV,,)- 



(B13) 
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